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ABSTRACT 

We report the detection of a strong Milky Way-type 2175 A extinction bump at 
z = 2.1166 in the quasar spectrum towards SDSS J121143.42-1-083349.7 from the 
Sloan Digital Sky Survey (SDSS) Data Release 10. We conduct follow up observations 
with the Echelle Spectrograph and Imager (ESI) onboard the Keck-II telescope and 
the Ultraviolet and Visual Echelle Spectrograph (UVES) on the VLT. This 2175 A 
absorber is remarkable in that we simultaneously detect neutral carbon (C i), neutral 
chlorine (Cl i), and carbon monoxide (CO). It also qualifies as a damped Lyman alpha 
system. The J1211-I-0833 absorber is found to be metal-rich and has a dust depletion 
pattern resembling that of the Milky Way disk clouds. We use the column densities of 
the C I fine structure states and the C ii/C i ratio (under the assumption of ionization 
equilibrium) to derive the temperature and volume density in the absorbing gas. A 
Cloudy photoionization model is constructed, which utilizes additional atoms/ions to 
constrain the physical conditions. The inferred physical conditions are consistent with 
a canonical cold (T ^ 100 K) neutral medium with a high density (n(H i) ^ 100 
cm“^) and a slightly higher pressure than the local interstellar medium. Given the 
simultaneous presence of C i, CO, and the 2175 A bump, combined with the high 
metallicity, high dust depletion level and overall low ionization state of the gas, the 
absorber towards J1211-1-0833 supports the scenario that the presence of the bump 
requires an evolved stellar population. 
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(LMC) exhibits this feature as well, but the canonical Small 
Magellanic Cloud (SMC) extinction curve does not show a 
bump. The trend of decreasing bump strengths in the MW, 
LMC, and S MC corresponds to progressively lower metal 
abundance s llFit zpatriclj |2004|l or d i fferen t radiative envi- 
ronments llCordon. Calzetti fc Witl] Il997l : iMattsson et al.l 

l2008ll . 


1 INTRODUCTION 

The most prominent feature in the interstellar extinction 
curve of the Milky Way (MW) is the broad absorption bump 
centered at rest-frame 2175 A. The Large Magellanic Cloud 
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Measuring the extinction curves in the Local Group 
by comparing stellar spectra is not feasible at larger dis¬ 
tances. Instead, more luminous background sources are re¬ 
quired to probe the gas and dust content beyond the local 
Universe. Quasars, one of the brightest populations in the 
Universe, have been used as background sources to study 
absorption line systems. The intervening quasar absorption 
lines serve as a powerful tool for studying physical condi¬ 
tions, chemical abundances and gas kinematics in interven¬ 
ing absorbers, such as damped Lyman-a absorbers (DLAs), 
Mg II absorbers, Ca ll and C iv absorbers. To date, a few 
dozens of the 2175 A bump detections have bee n reported in 


the s p ectra of intervening Mg ll absorber s (e.g. , lWan£_et_^ 


2004 


201 c 

; 

2003 



Srimiaiidet alJ 20081 : Jiang et al.l 12011 


Zhou et al 


Zhnn^_et_aL 2015[) and DLAs (e.g.. IWucknitz et al 


Junkkarinen et al. l2004l : IWang et al.l 2 OI 2 I 1 towards 


quasars. Most of the det ections are r eporte d in individual 
absorption line systems. Ijiang et al.l (l201lli systematically 
selected thirty-nine 2175 A absorbers at z ^ 1-1.8 asso¬ 
ciated with strong Mg ll absorption lines f rom the Sloan 
Digital Sky Survey (SPSS: [York et ^l2000ll Data Release 
3, which is the largest sample to date. Most bumps in this 
sample resemble the relatively weak 2175 A bump in LMC 
supershell rather than the strong bump observed in the 
Milky Way. The 2175 A bump featu re is also found in rare 
dust-rich and m etal-strong DLAs dProchaska et al] l2009l : 
IWang et al.ll20l3l . although DLAs are generally metal-poor 
system s with only a few perce n tage of the solar me tal abun¬ 
dance (IProchaska et al] I2OO3I : iKaplan et al.l I2OI0I I . Recent 
years have seen the rapid growth of studies of Gamma-ray 
burst (GRB) afterglows which are promising in probing even 
higher redshift absorbing systems. A few cases of GRB af¬ 
terglow spectra also reveal the 2175 A bump from inter- 
vening systems or from gas in the GRB host galaxies (e.g., 


Prochaska et al.l I 2 OO 9 I : lElfasdottir et al] I 2 OO 9 I : I Zafar et ahl 

2012f) . 

The nature of the 2175 A bump has been in- 
vestigated for almost 50 years since its first discovery 
(IStecheij Il965^ . iNoll et ^ (l2007ll found a correlation be¬ 
tween heavy reddeni ng and the presence of the bump. 
Severa l studies (e.g., iNoterdaeme et aP I 2 OO 9 I : Ijiang et al.l 
l2010bl i found that a high dust depletion level (i.e., [Fe/Zn] 
~ -1.5) is required to produce a prominent 2175 A bump. 
Simultaneous detections of metal lines, neu tral carbon 
llElfasdottir et al.l2009l: In oterdaeme_et_alj2009jL_and molec¬ 
ular gas ( Prochaska et al] 20091 : Noterdaeme et aT1 l2009l l 
with the 2175 A bump tend to suggest a correlation between 
those ingredients and the presence of the bump. The exact 
origin or carrier of the 2175 A bump remains an open ques¬ 
tion, although candidates have been suggested to be some 
form of graphitic car bon, most li kely polycyclic aromatic 
hydrocarbons ('PAHs: IPrainelliooiil L 

We initiated a project to search for 2175 A bumps in 
SDSS quasar spectra, yielding a large sample of 2175 A ab¬ 
sorbers (Zhang et al. in prep). From this survey, an interven¬ 
ing system towards SDSS J121143.42-1-083349.7 (hereafter 
J1211-|-083(f]) stands out as a remarkable 2175 A absorber 
with a significant absorption feature and rich metal lines. 


Along with the bump in the same absorber, we have simul¬ 
taneously detected G i fine structure lines, GO absorption 
bands, and molecular hydrogen hand in hand with neutral 
chlorine. 

G I fine structure transitions have been recognized 
and used to probe cold gas in high-redshift universe. The 
ground state of G l is split into three fine structure levels 
referred to as G l, G l*, and G l**. They are sensitive 
probes due to the fact that the energy separations between 
excited states and the ground state are only 23.6 K and 
62.4 K, respectively. Physical conditions of the absorber 
(i.e., temperature, density, and pressure) can be implied 
by the relative excitatio n of the G l fine structure states 
djenkins fc Shave] Il979fl . In the Milky Way, surveys of 
inter s tellar G l f ine structure excitations (Ijenkins fc Tripnl 
l200ll . l2007l . I 2 OIIII find that in the local interstellar medium 
(ISM) the gas pressure is mostly between 3 < log(P/k) <4 
cm“^ K. G I analysis has revealed the physical conditions 


in several high - z DL A s fe.g.. iG e. Bechtold fc BlackI Il997l : 
ISrianand et al.l l2005l : INoterdaeme et al.l 2001^) together 


systems fe.ff.. ISrianand et a 

J I 2 OO 8 I: IProchaska et alJ l2009l: 

Noterdaeme et al. 20091: 

Jorgenson. Wolfe & Prochaskal 

20in| ). Analvzing the excitation/de-excitation of the 


molecular and atomic species together provides detailed 
physical conditions of the absorbing gas and much needed 
insight into the properties of the intervening galaxy (i.e., 
its chemical abundances, dust content, and star formation 
rate). 

This paper examines the properties of this particular 
2175 A absorber, and is organized as follows. In Section O 
we describe target selection, observations, and data reduc¬ 
tion conducted on J1211-I-0833. In Section |3l we derive the 
bump strength, bump width and the extinction curve. In 
Section |4l we show the velocity profiles of absorption lines 
and column density measurements. Gas-phase metal abun¬ 
dances and dust depletion pattern are derived in Section O 
In the next section, we report the detections of G i, G1 i, 
and GO. We derive the physical conditions through the G l 
method and also construct a Cloudy photoionization model. 
The results are discussed and summarized in Section 0 


2 TARGET SELECTION, OBSERVATIONS, 
AND DATA REDUCTION 

J1211-I-0833 was originally selected from the tenth data re¬ 
lease (DRIO) of the S DSS-III’s Baryo n Os cillation Spectro¬ 
scopi c Survey (BOSS; Dawson et al.ll2ni3l : lEisenstein et al.l 
I2OIII: IRoss et al.ll2012lL It was targeted for spectroscopic 
observations on 2012 January 19. The BOSS spectra 
ijSmee et al.ll2013h cover a wavelength range of 3650-10400 
A with a spectral r esolution of R ~ 2000. The quasar red- 
shift is Zqso=‘2A83 (jBolton et al .11201^ : I Paris et al.ll2014l L A 
significantly broad absorption feature is clearly visible in the 
spectrum (Fig. [T]), which is later confirmed as the 2175 A 
extinction bump. The point-spread-function ugriz magni¬ 
tudes for J1211-b0833 are u = 21.10 ± 0.10, g = 19.63 ± 


^ J1211-1-0833 is used to name the quasar. We refer to 
the absorber as “the (intervening) absorber/system towards 


J1211-I-0833” or “the J1211+0833 absorber” throughout the pa¬ 
per. 
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0.01, r = 19.47 ± 0.02, i = 19.15 ± 0.02, and 2 = 18.41 ± 
0.03. 

We conducted follow-up observations of J1211-I-0833 
with the Echelle Spectrographs and Imager (ESI; 
ISheinis et al.l I 2 OO 2 II onboard the Keck II telescope in 
order to reveal at higher spectral resolution the detailed 
properties of the ISM towards the quasar J1211-I-0833. The 
spectra were obtained on 2013 March 8 with two exposures 
of 1800 s and 1500 s, using the 0.75” slit with a corre¬ 
sponding resolution of R ~ 5400 . The seeing conditions 
during observation varied from 1.5” to 2”. The data were 
reduced, calibrated, and combined by the ESI Echelle Data 
Reduction package Q The hnal combined spectrum covers a 
wavelength range of 3900 A to 11715 A with a median S/N 
~ 9 per pixel around 6500 A. The standard star G19B2B 
was used for relative flux calibration. 

.11211-1-0833 was also observed with the Ultraviolet and 
Visible Echelle Spectrograph (UVES; Dekker et al. 2000) 
mounted on the VLT on March 25 and April 24, 2014 (Pro¬ 
gram ID: 093.A-0126) under good seeing conditions (FWHM 
< 0.7 ”). Four 4350 s exposures were taken simultaneously 
with the blue and red arms. We used a 0.9 "slit and a pixel 
binning of 2 x 2, leading to a resolving power of R ~ 54000 
and a mean S/N of ~ 6 . The data were reduced using the 
UVES pipeline (Ballester et al. 2000). The individual expo¬ 
sures were scaled and combined together, weighted by the 
inverse variance in each pixel. The cosmic ray impact resid¬ 
uals were rejected at the same time. The actual wavelengths 
covered in the reduced spectra are from 3288 A to 6648 A 
with two gaps between 4525 - 4621 A and 5600 - 5647 A. 


3 THE 2175 A EXTINCTION BUMP 

To extract the properties of the 217 5 A extinction bump 
towar ds J1211-I-0833, we adopt the iFitzpatrick fc Massal 
lll990li parametrization (EM parameterization) of the op- 
tical/UV extinction curve. The extinction curve in the rest 
frame of the absorber can be described as a combination of 
a linear component, characterizing the underling extinction, 
and a Drude profile, modeling the potential 2175 A bump. 
The extinction curve is parameterized as 

A{\) = Cl-\-C 2 X-\-czD{x,xo,'^), (1) 

where x = A“^, and the Drude profile D{x,xo,^) has the 
form of 

( 2 ) 

where xo and 7 are the bump peak position and bump 
width (FWHM) of the Drude profile, respectively. The 
UV linear component is set by the slope C 2 and inter¬ 
cept Cl. The derived extinction curve is not the absolute 
extinction curve normalized by E{B — V) but in a rela¬ 
tive sense. The strength of the bump can be defined using 
these parameters: Atump = 7103 /( 27 ) measures the area of 
the bump 0. We adopt the SDSS DR7 composite quasar 

^ http://www2.keck.hawaii.edu/inst/esi/ESIRedux/ 

® The A^„^„ defined in this p aper is slightly different from that 
in [Fitzpatrick &: Massal ll2007ll since their parameterized extinc¬ 
tion curves are normalized by E{B — V). Therefore = 


spectrum lijiang et al.ll201l]i as the intrinsic spectrum for 
J1211-I-0833, w hich is an update of the SDSS quasar com pos- 
ite spectrum of IVanden Berk. Richards fc Bau^ (1200 il l by 
median combining 105783 quasar s pectra in the SDSS DR7 
I Schneider. Richards fc Halil l2010l i. The composite spec¬ 
trum is reddened using the FM parameterized extinction 
curve at 2a6s=2.1166 to form the model spectrum. The ab¬ 
sorber’s redshift is based upon the strong absorption lines 
such as Mg II in the Keck spectrum. To optimize the contin¬ 
uum fitting, strong emission lines such as Lya, Si iv, C iv, 
C III, and Mg ll and known absorption lines are masked out. 

We compare the model spectrum with the ob¬ 
served BOSS spectrum which has been corrected 
for Galactic extinction by using the dust map of 
iSchlegel. Finkbeiner fc DavisI (1 1998 1 to find the best- 
fit. Fig. [T](a) demonstrates the fitting result with the black 
spectrum being the observed spectrum and the red curve 
representing the best-fit model. The best-fit parameters are 
Cl = -0.13, C 2 = 0.43, C 3 = 1.78, xo = 4.6, 7 = 1.47, and the 
bump strength Abump = 1.90 with = 1.21. In contrast, 
the green curve shows the composite quasar spectrum 
reddened by the best-fit linear component without an 
extinction bump, which explicitly illustrates the existence 
of a strong 2175 A bump. As the intrinsic spectrum of 
J1211-I-0833 is unknown, we should also consider the 
possibility of an intrinsically red or blue quasar spectrum 
that potentially provides a better fit. We select the bluest 
and reddest 20% quasars respectively from the SDSS DR7 
quasar sample according to spectral index. The blue and 
red quasar templates are constructed by median combining 
those selected quasar spectra. We perform the same fitting 
method using the blue and red templates (Fig[2]) as the 
intrinsic quasar spectra for J1211-I-0833, respectively. 
Comparing with the composite spectrum, an intrinsically 
blue or red quasar template provides an inferior fit in 
terms of representing the observed spectrum, which justifies 
the usage of the composite spectrum as a better model 
representative for J1211+0833. The bump strength and 
width increase to Abump = 2.42 and 7 = 1.67 in the case of 
the blue template and drop to Abump = 1.15 and 7 = 1.11 
for the red template. 

The strength of the 2175 A bump can appear to be 
stronger or weaker due to the intrinsic variation of quasar 
spectra, which in some cases imit ates an extinction bump 
dPitman. Clayton fc Gordonll200Clf l. We therefore utilize the 
simulation technique developed by Ijiang et al.l (l2010al l to 
set the criteria for the significance of the bump. The first 
step is to select a control sample of SDSS quasar spectra 
(sample size ~ 1000) at a similar redshift to the quasar of 
interest. For J1211-I-0833, we select quasars with redshifts 
in the range of 2qso-0.05< 2 < Zqso+0.05 and i-band S/N 
^ 5 as its control sample. Then we fit each of them by red¬ 
dening the composite quasar spectrum with a parameterized 
extinction curve at the redshift of the bump absorber of in¬ 
terest. The parameters xo and 7 in the Drude profile are 
held fixed to the best values derived from the fitting of the 
quasar of interest. The expected bump strength distribution 
is a Gaussian assuming random fluctuations of continuum 


E (B — V) X A*__ wh ere A^^^^ is the bump strength defined 

in lFitzpatrick fc M^sal ll2007l V 
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of each spectrum in the control sample. Bumps that signif¬ 
icantly deviate from the Gaussian distribution are defined 
to have statistical significance. The histogram in Fig. [T] (c) 
shows that the bump towards J1211-I-0833 lies at a statis¬ 
tical confidence level greater than 5a indicated by the red 
arrow. 

The best studied 2175 A bumps are those along lines 
of sight in our Milky Way and towards the LMC. We com¬ 
pare the bump strengths in the — 7 space (Fig. [3]) 

where the black filled circ les are the bumps in the MW 
dFitzpatrick fc Massa 2007ll and the g reen circles are LMC2 
supershell bumps 1 Gordon et al.llioosl l . The J1211+0833 ab¬ 
sorber possesses a typical MW-type bump strength but is 
wider than the majority. Does this indicate a mature envi¬ 
ronment resembling the MW has enabled the formation of 
such a strong extinction bump at the absorber’s redshift of 
a = 2.12? We resort to the column density measurements, 
metal abundances, and dust depletion pattern analysis to 
explore the physical and chemical environment in the ab¬ 
sorber. 


4 ABSORPTION LINES AND COLUMN 
DENSITY MEASUREMENTS 

4.1 Metal lines in Keck/ESI 

Numerous narrow absorption lines are detected in the 
Keck/ESI spectra of J1211-I-0833, enabling measurements 
of gas-phase column densities of many ions. The detected 
metal lines are listed in Table [T] and the absorption lines of 
interest are plotted in velocity space (i.e., normalized flux 
vs. relative velocity) in Fig. U) 

To measure column densities of metal absorption lines 
in the 2175 A absorption system, we normalize the ob¬ 
served Keck spectrum locally by the quasar continuum, 
which is fitted by polynomials with the absorption feature 
masked during the fitting. In the framework of column den¬ 
sity mea surements, the Apparent Optical Depth Method 
(AODM; ISavage fc Sembadil 1 199 il l is often applied to re- 
solved spectra and the Curve o f Growth (COG) method 
ll Jenkins. Savage fc Spitz'll 1986l l is commonly used in unre¬ 
solved cases. As the Keck/ESI spectrum for J1211-I-0833 is 
partially resolved, we adopt both methods to derive the col¬ 
umn densities as a consistency check. The AODM is based 
upon the connection between optical depth and column den¬ 
sity: the column densities are measured by integrating the 
optical depth over the velocity range of the spectra covering 
the detected transitions. The column densities derived from 
this method are lower limits for the saturated lines. For the 
COG method, the Fe ii AA 1608, 2344, 2374, 2382, 2586 and 
2600 absorption lines, arising from the ground level Ee"*", 
are used to construct the COG in Fig. (5] They all have the 
same column density A’(Fe ll) and the Doppler parameter 
b. The rest frame equivalent widths (EWs) are measured by 
integrating over a typical velocity interval of [-400, 400] km 
s“^ to include the absorption feature as well as to avoid con¬ 
tamination from other absorption. The best-fit COG is ob¬ 
tained by finding the best match between theoretical COG 
calculated with various b values and the measured EW of 
the six Fe ll absorption lines. We explore the b — N{Fe ll) 
parameter space based on minimization to find the best 


fit parameters log N[Fe ll) = 14.70 ± 0.07 cm“^ and b = 
47.5 ± 3.0 km s“^. The uncertainties quoted refer to the la 
confidence level. The column densities of other atoms/ions 
are derived based on the best-fit COG and the measured 
EWs. For species with multiple absorption lines, the col¬ 
umn densities are derived by simultaneously fitting all the 
available transitions. The Naodm and Ncoa values are in 
good agreement for the unsaturated lines. 


4.2 Metal lines in VLT/UVES and velocity 
structure 

The Doppler parameter derived from metal lines in the 
Keck/ESI spectrum implies there is blending, which is con¬ 
firmed by the high resolution (5.5 km s“^) VLT/UVES 
spectrum. The metal absorption lines are modeled by Voigt 
profiles using the VPFIT package (version 9.5) 0. VPFIT 
utilizes a x^-minimization algorithm for multi-component 
Voigt-profile fitting to derive absorption redshifts, Doppler 
parameters, and column densities taking into account the in¬ 
strumental resolution. A set of different singly ionized lines 
covered in VLT/UVES (i.e., Zn ll, Fe ll, Ni ll, Cr ll) are fitted 
simultaneously assuming they are kinematically associated 
with the same gas cloud. The redshift and b parameter for 
each component are tied to each other. We find that a mini¬ 
mum of nine velocity components are required to optimally 
fit the absorption features in strong Zn ll and Fe ll lines (Fig. 
O. We number the components from blue to red as compo¬ 
nents 1 to 9 at the following velocities relative to Zaba = 
2.1166 for u = 0 km s“^(the center identified in Keck/ESI): 
V -52 km s“^, -31 km s“^, -19 km s“^, -8 km s“^, -|-5 
km s“^, -1-26 km s“^, -1-38 km s“^, -1-84 km s“^, and +99 
km s“^. Ni II and Cr ll line profiles are attributed to the pri¬ 
mary components 3 and 4 at 2 = 2.116404 and 2 = 2.116520, 
respectively. The nine components span a velocity interval 
of Av ~ 150 km s“^. Note that the correspondence be¬ 
tween Av and b (from the cu rve of growth on ESI) ma tches 
very well the relation from iNoterdaeme et ^ (l2014l ~): Av 
= 2.216e// + 0.02b^ff. The profiles exhibit an edge-leading 
asymmetry with the strongest feature in the blue end. The 
asymmetric shape is consiste nt with the prediction by a sim¬ 
ple model of a rotating disk JProchaska fc Wolfelll997l l. 

The column densities in each component and associated 
errors are provided by VPFIT and are reported in Table [S] 
The total column densities for each ion summed over all the 
components, log A(Zn ii) = 13.58 ± 0.06, log N(Fe ii) = 
14.77 ± 0.03, log A(Ni ii) = 13.62 ± 0.08, and log iV(Cr ii) 
= 12.78 ± 0.14 cm“^, are consistent with the measurements 
from Keck/ESI. Potential contamination to Zn ll column 
density from Cr ll A2026 and Cr ll A2062 lines is taken into 
account in the simultaneous fitting. Note that the b values 
obtained from the fit can be smaller than the velocity res¬ 
olution of the instrument because the Doppler parameters 
are constrained by simultaneously fitting several transitions 
with a range of oscillator strengths (i.e., probing different 
regimes along the curve of growth). The real resolution can 
be better than the nominal resolution we used in the fitting 
because the seeing is less than the slit width. This may well 
explain the small 6-values we obtained. 


^ http://www.ast.cam.ac.uk/~rfc/vpfit.html 
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Figure 1. The 2175 A bump fitting results for the absorber towards J1211+0833. (a) The black spectrum is the observed quasar 
spectrum retrieved from BOSS DRIO. The red curve is the best-fit model with the extinction bump whose central position is indicated 
by the red arrow and the blue arrow points to the Mg ii absorption lines. The green curve is the best-fit model reddened by the linear 
component only, (b) The best-fit extinction curve with the constraint that the extinction should be zero at A —>• oo. (c) The bump 
strength distribution of the control sample fit by a Gaussian curve in blue. The vertical dashed blue line indicates 5cr and the red arrow 
denotes the bump strength towards J1211-h0833. 
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Figure 2. The 2175 A bump fitting results using the blue (left) and red (right) quasar templates. The black spectrum is the observed 
quasar spectrum retrieved from BOSS DRIO. The red curve is the best-fit model with the extinction bump the central position of which 
is indicated by the red arrow and the blue arrow points to the Mg II absorption lines. The green curve is the best-fit model reddened by 
the linear component only. The 2175 A bump is present even in the extreme case of the red template. 
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Figu re 3. The co mparison of 2175 A b ump strengths Ahump bump widths 7 with extinction bumps in the MW l|FitzDatrick Massal 
l2007t) and LMC2 llGordon et alj[2003h supershell denoted with black and green circles, respectively. The absorber towards J1211+0833 
is denoted as the red star. The error bars reflect the uncertainty induced by the blue and red template fitting. 


Table 1. Column density measurements of the system towards 31211+0833 using the Keck/ESI data. The central wavelengths of the 
transitions a nd the equivale nt widths are reported in the absorber rest frame. We adopt the oscillator strengths, /, from the Atomic Data 
compiled bv iMortonI l|20Q3l l. Saturated lines have lower limits placed on their column densities derived from the AODM. Undetected 
lines have 3cr upper limits on equivalent widths and hence column densities from the COG method. 


Ion 

Transition 

(A) 

/ 

EW 

(A) 

log Naodm 
(cm“^) 

log Ncog 

(cm-2) 

Fell 1608 

1608.4511 

0.0577 

0.52 ± 0.13 

14.74 ± 0.07 

14.70 ± 0.07 

Fell 2344 

2344.2139 

0.1140 

1.01 ± 0.13 

14.52 ± 0.03 


Fell 2374 

2374.4612 

0.0313 

0.50 ± 0.06 

14.61 ± 0.06 


Fell 2382 

2382.7652 

0.3200 

1.36 ± 0.12 

14.59 ± 0.01 


Fell 2586 

2586.6500 

0.0691 

1.04 ± 0.06 

14.62 ± 0.03 


Fell 2600 

2600.1729 

0.2390 

1.33 ± 0.10 

14.73 ± 0.01 


Mgl 2852 

2852.9631 

1.8300 

0.88 ± 0.10 

12.98 ± 0.03 

13.08 ± 0.10 

Mgll 2796 

2796.3553 

0.6155 

2.43 ± 0.06 

>14.20 

15.97 ± 0.09 

Mgll 2803 

2803.5324 

0.3058 

2.10 ± 0.05 

>14.47 


Sill 1526 

1526.7070 

0.1330 

1.05 ± 0.06 

>14.95 

15.65 ± 0.07 

Sill 1808 

1808.0129 

0.0021 

0.20 ± 0.04 

15.64 ± 0.10 


Nill 1751 

1751.9157 

0.0277 

0.03 ± 0.01 

13.68 ± 0.69 

13.61 ± 0.18 

Znll 2026 

2026.1370 

0.5010 

0.37 ± 0.07 

13.43 ± 0.08 

13.51 ±0.09 

Znll 2062 

2062.6604 

0.2460 

0.32 ± 0.07 

13.61 ± 0.11 


Mnll 2576 

2576.8770 

0.3610 

0.27 ± 0.05 

13.14 ± 0.11 

13.09 ± 0.06 

Mnll 2594 

2594.4990 

0.2800 

0.20 ± 0.04 

13.10 ± 0.14 


Mnll 2606 

2606.4620 

0.1980 

0.12 ± 0.04 

12.69 ± 0.48 


CrII 2056 

2056.2569 

0.1030 

<0.25 


<13.95 

CrII 2062 

2062.2361 

0.0759 

<0.43 



CrII 2066 

2066.1640 

0.0512 

<0.28 



Aim 1854 

1854.7184 

0.5590 

0.45 ± 0.10 

13.76 ± 0.04 

13.72 ± 0.07 

Aim 1862 

1862.7910 

0.2780 

0.34 ± 0.04 

13.79 ± 0.08 


SilV 1393 

1393.7602 

0.5130 


>14.49 


SilV 1402 

1402.7729 

0.2540 


>14.71 


CIV 1548 

1548.1950 

0.1908 


>15.02 


CIV 1550 

1550.7700 

0.0952 


>15.27 
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Figure 4. Metal absorption lines detected in the Keck/ESI for the absorber at 2:=2.1166 towards J1211-I-0833. 
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log (Na) 


Figure 5. The COG using six Fe ii absorption lines from Keck/ESI. The three curves represent the best-fit and corresponding la 
confidence level. The column densities of other atoms/ions are derived based on the best-fit COG and the measured EWs. 
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Figure 6. Simultaneous multi-component Voigt-profile fitting 
with VPFIT on the VLT/UVES data. The original spectra are 
in black and the fitted profiles are in red. A minimum of 9 ve¬ 
locity components indicated by the vertical lines are required to 
optimally fit the absorption features. 


4.3 Lyman-a profile fitting 

The Lya absorption line of the system towards J1211-1-0833, 
covered by the VLT/UVES spectrum, is fitted by a Voigt 
profile using xjitdla from the XIDL package 0. We adopt 
the best-fit Voigt profile (Fig. ^ with a neutral hydrogen 
column density of log A'’(H l) = 21.00 ± 0.20 cm“^. The 
error reflects the systematic uncertainty induced by the con¬ 
tinuum level plus small uncertainty in the fitting. We also 
utilize VPFIT to perform the Ly-a profile fitting which is 
based on normalized flux. We use both as a sanity check. 
The resulting best-fit log V(H l) is 0.05 dex higher and well 
within the estimated error. The sys tem towards J121 1-h0833 
is therefore classified into a DLA (IWolfe et al.lll98^ . With 
the hydrogen column density available, we are able to derive 
the absolute metal abundances. 


5 GAS-PHASE ABUNDANCES AND DUST 
DEPLETION PATTERN 

The column densities of ions measured from absorption lines 
refer to the species in the gas phase. We assume the ioniza¬ 
tion correction of metals is negligible due to the shielding 
effect of large column density of neutral hydrogen from ion¬ 
izing photons. The total elemental column densities can be 
derived from the dominant ionization state, i.e., the lowest 
energy species where the ionization potential is above 13.6 
eV. The gas-phase abundances reveal the underlying nu- 
cleosynthetic abundances modified by differential depletion 
pattern of elements. Heavy element abundances are mea¬ 
sured relative to solar values, as [X/H] = log N{X)/N{H) 


® http://www.ucolick.org/~xavier/IDL/index.html 































































J1211 9 



3740 3760 3780 3800 3820 3840 

Wavelength (A) 


Figure 7. DLA profile fitting on the VLT/UVES spectrum (grey). The spectrum plotted in black is the original spectrum smoothed 
using 5-pixel boxcar and the red dotted line represents the associated flux errors. The blue line is the best-fit Voigt profile with a hydrogen 
column density of log Af(H l) = 21.00 cm“^ and the estimated 1 cr uncertainty is denoted by the blue shade. The blue dashed line is the 
associated continuum. 


- log {X/H)q assuming that Af(H) = N{R l). The resul¬ 
tant metal abundances are listed in Table [S] The Zn ll col¬ 
umn density indicates that this is a metal-strong DLA (lo g 
Af(Zn ll) (S 13.15) as defined bv iHerbert-Fort et al.l (l200fil ). 

FigEldemonstrates metallicity measured in [Zn/H] ver¬ 
sus hydrogen column density log Af(H l) in different ab¬ 
sorbers. The crosses are the DLAs with [Zn/H] mea¬ 
surements in IProchaska et al.l ||2007|) and the ope n cir¬ 
cles are clouds in the MW by iRoth fc Blades! (Il995l ) . The 
J1211-I-0833 absorber, with a moderate hydrogen column 
density (in the DLA regime), distinguishes itself by having 
a metallicity ([Zn/H] = -0.07 ± 0.21) comparable to that of 
the MW clouds and above all the DLAs in this sample. 

In Fig. [9] we show the dust depletion pattern of the 
absorber towards J1211-1-0833 compared to that of the 
MW cool and warm diffu se disk clouds toward Ophiuchi 
dSavage fc Sembachlll99(3 ). A very high relative abundance 
of [Zn/Fe] = 1.74 indicates the metals are heavily depleted 
onto dust grains. The dust depletion level of the absorber to¬ 
wards J1211-I-0833 lies between that of the warm disk clouds 
and the cold disk clouds with an enhanced Zn abundance. 
The high depletion level supports the existence of a large 
number of dust grains and thus the observed strong 2175 A 
dust extinction bump. 

6 PHYSICAL CONDITIONS THROUGH 

ATOMIC AND MOLECULAR LINES 

6.1 Cl profiles and excitation temperatnres 

The strongest multiplets at 1560 A and 1656 A are detected 
in the Keck/ESI spectra, and many close fine structure tran¬ 
sitions are heavily blended. We use the VPFIT package to 
simultaneously fit as many C l multiplets as possible in the 
UVES data: the strongest multiplets at 1560 A and 1656 


A that fall redward of the Lya forest and the multiplet at 
1328 A. The other multiplets at bluer wavelengths are re¬ 
jected due to low signal-to-noise ratio and contamination. 
Fig. [To] demonstrates the three multiplets in black, overlaid 
with our fit in red, and the C l fine structure transitions are 
marked with the vertical lines. Each fine structure state is 
denoted by a different color (C l: blue, C T: orange, C T*: 
green). As shown in the figure, nine velocity components are 
required to produce a satisfactory fit, which is similar to the 
velocity structure of low ionization lines obtained from Zn ll, 
Fe II, Ni II, and Cr ll. The components are located at u ~ 
-51 km s“^, -28 km s“^, -17 km s“^, -6 km s“^, -1-5 km 
s“^,-1-20 km s”*^,-1-39 km s“^,-1-79 km s“^, and-1-99 km 
s“^. A multi-component fit on a neutral chlorine line (Fig. 
dH), Cl I A1347, reveals a similar velocity structure with C l, 
which is ex pected due to their similar ionization potentials 
l|Jud[l97j) (more discussion on chlorine in Section (6.411 . 

The resultant C l, C T, C T* column densities in each 
component are summarized in Table |4| We also derive the 
excitation temperature (Tex) of C l according to the Boltz¬ 
mann equation. 


Ne/Ng = ge/ggexp{-AEeg/kTex), (3) 


where and Ng are column densities in the excited and 
ground state levels, and g^ and gg are the corresponding 
statistical weights, k is the Boltzmann constant. AE^g (if 
given in K) is 23.6 K for the J—1 to 0 transition in C l. 
The excitation temperatures derived from the C l and C T 
states in each component are all consistent with and higher 
than the cosmic microwave background (CMB) temperature 
Tomb = 8.5 K at the absorber’s redshift. T = Tq(1 -\- z) 
where To = 2.725 K llMather et al.l[l999l ). 
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Figure 8. Metallicity comparison with MW clouds llRoth fc Bladeslll99Rh and normal DLAs llProchaska et aLll2007h in terms of [Zn/H]. 
The J1211+0833 absorber, with a moderate hydrogen column density, distinguishes itself by having a metallicity ([Zn/H] = -0.07 ib 0.21) 
comparable to that of the MW clouds and above all the DLAs in this sample. 



Figure 9. Dust de pletion pattern of the abs orber towards J1211-1-0833 in comparison with that of the warm and cool diffuse disk clouds 
toward Ophiuchi llSavage Sembachlll99^ . A very high relative abundance of [Zn/Fe] = 1.74 indicates the metals are heavily depleted 
onto dust grains. The dust depletion level towards J1211-I-0833 lies between that of the warm disk clouds and the cold disk clouds with 
an enhanced Zn abundance. 
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Figure 10. C I velocity structure (nine components) derived from simultaneous fitting of three multiplets: C I 1560, C I 1656, and C I 
1328 in the VLT/UVES data. Black is the data and the overall fit is in red. The different fine structure transitions are color-coded by 
the vertical lines-C I: blue (dashed lines), C I*: orange (solid lines), C I**: green (dotted lines) 


6.2 Deriving physical conditions through C l fine 
structure lines and the C ll/C l ratio 

In order to derive the physical conditions in 
the absorbing gas, we e mploy the method by 
Ijorgenson. Wolfe fc Prochaskal ll2010l ) which simultane¬ 
ously constrains both the volume density and temperature 
of the absorbing gas without making assumptions of tem¬ 
peratures. We invoke the assumptions of steady state and 
ionization equilibrium. Various excitation and de-excitation 
mechanisms are taken into account: spontaneous radiative 
decay, direct excitation by the CMB, UV pumping due to a 
radiation field, and collisional excitation and de-excitation. 
The volume densities and temperatures are constrained 
through two steps: first solving the steady state equation 
by calculating the ratio of each excited state relative to 
the ground state, and then further constraining the C l 
solutions with the ionization equilibrium invoked by the 
C ll/C I ratio. 

The rate of populating and de-population of state i is 
given by 

'’T'j + BjiUji + Fji + J2k = /-N 

+ BijUij + Fij + J2k <1%) 


where Aij is the probability of spontaneous decay , Bij is 
the probability of stimulated emission, Uij is the energy den¬ 
sity of the cosmic microwave background radiation field, Fij 
is the indirect excitat ion rate due to fluorescence defined by 
ISilva fc ViegasI ll20nil ). In addition, the excitation and de¬ 
excitation terms due to collisions are included. The quantity 
n* is the volume density of the collision partner k, where k 
can be atomic hydrogen, electron, and proton (densities de¬ 
noted as n(H l), rie, and Up) in the case of a DLA and qfj 
= < (TU > is the collision rate coefficient, where a is the 
cross section and v is the relative speed between the col¬ 
lision partners. The reverse rates are calculated using the 
assumption of detailed balance. Taking all of these excita¬ 
tion and de-excitation mechanisms into account, the ratios 
of the upper to lower fine structure level populations of C° 
become a function of redshift, temperature, neutral hydro¬ 
gen volume density, and electron density. We first consider 
the UV pumping due to a radiation field that has a typical 
intensity of Jk ~ er gs cm~® s~^ Hz~^ sr“^measured 

in damped Ly-a systems I Wolfe et al.ll2004h and an elec¬ 
tron density of rie = 10“® n(H l) w hich is consistent with 
the typical values found in DLAs (|Srianand et al.l l2005l : 
iNeeleman, Prochaska fc WolfejEoiRh . For now, we use typi¬ 
cal values but note that they could actually differ strongly 
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Table 2. Doppler parameters and column densities for each com¬ 
ponent measured using the VLT/UVES data. 


Comp. Zabs lofi ^ (km s log N (cm 


1 

2.116061 Zn II 

5.58 

12.09+0.07 


Fe II 

10.60 

13.79+0.04 


Ni II 

10.60 

12.88+0.15 


Cr II 




2 

2.116277 Zn II 

8 .20+1.12 

12.49+0.09 


Fe II 

8.87 

13.93+0.07 


Ni II 

8.66 

12.67+0.28 


Cr II 




3 

2.116404 

Zn II 

Fe II 

Ni II 

Cr II 

8.38+1.28 

9.06 

8.84 

9.39 

13.07+0.05 

14.25+0.08 

13.32+0.08 

12.53+0.13 

4 

2.116520 

Zn II 

1.80+0.68 

12.91+0.23 



Fe II 

1.95 

13.36+0.30 



Ni II 

1.90 

11.79+1.37 



Cr II 

2.02 

12.27+0.17 

5 

2.116648 

Zn II 

12.02+1.18 

12.81+0.03 



Fe II 

13.01 

13.86+0.04 



Ni II 

12.69 

11.70+2.58 



Cr II 

13.48 

11.47+1.64 

6 

2.116875 

Zn II 

3.74+0.71 

12.04+0.10 



Fe II 

4.05 

13.74+0.06 



Ni II 

3.95 

12.61+0.21 



Cr II 




7 

2.116998 

Zn II 

9.53+1.33 

12.35+0.06 



Fe II 

10.31 

13.61+0.06 



Ni II 

10.06 

12.19+0.75 



Cr II 

10.69 

11.68+0.99 

8 

2.117469 

Zn II 

Fe II 

18.94+1.93 

20.49 

12.46+0.06 

13.78+0.05 


Ni II 
Cr II 


9 

2.117628 Zn II 

1.62+1.53 

11.97+0.12 


Fe II 

1.75 

12.98+0.14 


Ni II 

1.71 

12.30+0.36 


Cr II 




total 

Zn II 

13.58+0.06 


Fe II 

14.77+0.03 


Ni II 

13.62+0.08 


Cr II 

12.78+0.14 


Table 3. Gas-phase abundances of the absorber towards 
J1211-I-0833 using Keck/ESI and VLT/UVES. 



J1211+0833 

source 

Solar abundances 

[Zn/H] 

-0.07+0.21 

VLT/UVES 

-7.35 

[Si/H] 

-0.98+0.23 

Keck/ESI 

-4.45 

[Fe/H] 

-1.81+0.21 

Keck/ESI 

-4.49 

[Mn/H] 

-1.39+0.23 

Keck/ESI 

-6.47 

[Ni/H] 

-1.63+0.22 

VLT/UVES 

-5.75 

[Mg/H] 

> -2.11 

Keck/ESI 

-4.42 

[Cr/H] 

-1.90+0.24 

VLT/UVES 

-6.32 



Relative Velocity (km s ’) 


Figure 11. A multi-component fit on Cl i A1347 from 
VLT/UVES. The resultant total neutral chlorine column density 
is log N{Cl i) = 13.95 zb 0.10 cm-^. 


in this peculiar system. We discuss the sensitivity to the 
radiation field and electron fraction afterwards. 

We find the solution in the (/I, /2) plane where /I = 
n(C l*)/n(C l)to£; /2 = n(C i**)/n(C l)£o£5 3jlici Tl(^C l)£ot 
= n(C l) + n(C I*) +n(C l**). The values of (/I, /2) for 
each component are given in Tabled We show an example 
case of component 3 of the absorber towards J1211+0833 in 
Fig. 1121 According to the Zn ll column density distribution 
in each component, component 3 shows the highest column 
density with about 30% of the total. We refer to the Zn ll 
column density not only because its velocity profile is the 
best measured one available to us but also the fact that Zn 
is non-refractory and non-depleted such that it is likely to 
trace the underlying relative distribution of H l or other low 
ions in each cloud. We generate theoretical tracks: one for 
each temperature from T = 10 to 10^ K in steps of 0.1 dex; 
for each track n(H l) ranges from 10“®'® to lO'^'^ cm“® in 
steps of 0.02 dex. The orange/blue polygon, which is deter¬ 
mined by the 1 a/2 a error bars for the column densities, 
encloses the region where 1 cr/2 a solutions are accepted 
(Top Left). n(H l) is well constrained for each temperature, 
but temperature is still a free parameter based on C l fine 
structure lines only (Top Right). 

We further constrain the C l solutions by introduc¬ 
ing the C ll/C I ratio under the assumption of ionization 
equilibrium. The ionization equilibrium in this case can be 
written as rie n{C ll) a = n(C l) F, where a is the re¬ 
combination coefficient of element C'*' to C° and F repre¬ 
sents the ionization rate which is proportional to the ra¬ 
diation field intensity. The recombination coefficient a in- 
cludes the effects of radiative plus dielectric recombination 
dShull &L van Steenbers 1982ll and grain-ass isted collisional 
recombination ( Weingartner fc Drainel[200lf) . We cannot di¬ 
rectly measure the N{C ll) because the resonance lines 
of C II are too saturated. We instead me a sure N(C ii) 
by proxy using N{Si ll) as in I Wolfe etlil] ll2004l : [C/H] 
= [Si/H] + [Fe/Si] int where the intrinsic (nucleosynthetic) 
ratio [Fe/Si]int = —0.2 for a minimal depletion model or 
[Fe/Si]int = 0.0 for a maximal depletion model. We adopt 
an average [Fe/Si] int ratio of the two depletion models and 
account for the errors induced by choice of depletion mod- 
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els. We use A^(H l) and N{C ll) as global measurements 
and apply them to individual C l-bearing clouds. Since Ly- 
a is damped, the velocity structure associated with each H l 
cloud is unknown. For C ll, we cannot measure an accurate 
column density in each velocity component unless an unsat¬ 
urated C II or Si II line is covered but unfortunately Si ll 
1808 lies in the gap of the VLT/UVES spectrum. We do not 
attempt to perform strict individual component analysis in 
this regard. 

We plot log(C ll/C l) as a function of n(H l) to demon¬ 
strate the constraints imposed by ionization equilibrium 
(shown in the lower left panel of Fig. I12II . The measured log 
N{C ll/C l) = 2.46 ± 0.15 for component 3 of C l-bearing 
clouds is indicated by the horizontal dashed lines. The final 
solutions are re-plotted in the temperature vs. n(H l) dia¬ 
gram (Bottom Right panel). In the case of component 3, the 
density is constrained in the range 105 < n(H l) < 166 cm“® 
at the 1 a level and 91 < n(H l) < 191 cm“® at the 2 a level; 
the temperature is constrained to be 63 K < T < 200 K (1 
a) and 40 K < T < 316 K (2 a). The constrained pressure 
is 3.66 < log (-P/k) < 4.66 cm“® K (2 a). We can compare 
our results to those found in the local ISM. Jenkins & Tripp 
(2011) find by using C l fine structure lines that the cold 
neutral medium in the local ISM has an average pressure of 
log (P/k) = 3.58 ± 0.18 cm"® K. 

We further estimate the characteristic size and mass of 
the absorbing cloud. The size can be estimated by the ratio 
of N(H l) to n(H l): I = N(H l)/n(H l) ~ 2 - 4 pc. This is es¬ 
sentially an upper limit in a cold component because A^(H l) 
is distributed over all components. Therefore C l traces very 
dense pockets of very cold gas at slightly higher pressure. 
The associated gas cloud has a mass M = mpriHil^ (rtip 
is proton mass) 3110 - 9000 solar masses. A rough esti¬ 
mate of the mass of the galaxy is Mgai = u^R/G 5 x 10®° 
( lokpc ) ^0’ "'here v is the velocity span of the components, 
R is the radius, and G is the gravitational constant. 

The above results are based on the typical values of the 
radiation held intensity Ji, and electron fraction in DBAs. At 
those high hydrogen volume densities, the relative popula¬ 
tions of G I hne structure states are dominated by collisional 
excitation and de-excitation and the UV pumping has negli¬ 
gible effect. Whereas the G ll/C l ratio can be largely driven 
by the UV pumping rate as well as the electron fraction 
as G ll/C I is directly proportional to F/ue. We demon¬ 
strate how the theoretical C ll/C l ratio moves in the log 
N{C ll/C l) vs. n(H l) plane (Fig. I13II with a grid of Jv 
from the strength of the Haardt -Madau background (2.54 
X 10~^°: iHaardt fc Madaul [199611 to 10" ergs cm ^ s ® 
Hz"® sr"®indicated by different colors. The electron frac¬ 
tion is hxed to rie = 10"® n(H l). All the tracks are the 2 a 
results. Again, the allowed solutions are constrained by the 
measured C ll/C l ratio. The case of = 10"®® is ruled 
out by the measured value. The allowed solutions are con¬ 
strained in the range 40 - 5200 cm"®. Fig. 1141 shows the 
variation due to a grid of possible electron fractions (10"® 
to 10"®) with a fixed Ju of 10"®° ergs cm"° s"® Hz"® sr"®. 
The range of n(HI) is restricted to be 10 - 360 cm"®. Only 
by having good knowledge of these two quantities can we 
tightly constrain the temperature and density. 

The other uncertainty is induced by applying the total 
N{C ll) to individual G l-bearing clouds. Only a part of the 
total N{C ii) is actually associated with each C i cloud. The 



n(Hl) (cm‘^) 


Figure 13. The sensitivity of C Il/C I ratio to the radiation 
field intensity in the example case of component 3. Theoretical 
tracks corresponding to each J^, are marked by different colors. 
The electron fraction is held fixed. 



Figure 14. The sensitivity of C Il/C l ratio to the electron frac¬ 
tion in the example case of component 3. Theoretical tracks corre¬ 
sponding to each are marked by different colors. The radiation 
field intensity is held fixed. 

adopted G ii/C i ratio serves as an upper limit. A decrease in 
the amount of N{C ll) associated with the G l cloud would 
result in a higher volume density and lower temperature and 
therefore an approximately constant pressure. 

6.3 Cloudy photoionization modeling 

The physical conditions derived from the above method are 
average properties of the absorbing cloud or the same condi¬ 
tion throughout the cloud. This approach is sensitive to the 
adopted electron density and radiation field intensity which 
can be poorly const rained. We further construct a Cloudy 
llFerland et al.llioisl : version 13.03) photoionization model of 
this absorber by which we can utilize observed quantities in 
addition to C l fine structure lines. We model the absorber 
as plane-parallel slabs of gas of constant density exposed to 
a radiation field. We consider the spectral energy distribu¬ 
tion of the Milky Way ISM, table ISM built in Cloudy. The 
intensity of the ISM radiation can be changed by a scale 
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Figure 12. Top Left: fl vs. /2 for component 3 of the J1211+0833 absorber. The measured value is denoted by a red asterisk and the 
1 (7 (2 (t) error polygon is in orange (blue). The black curves are the theoretical tracks which cover possible ranges of temperatures and 
hydrogen volume densities. Top Right: Temperature vs. n(H i) for the allowed solutions (1 a in orange and 2 cr in blue) constrained by 
the C I fine structure lines only. Bottom Left: log A^(C Ii/C i) vs. n(H i) for the allowed solutions further constrained by introducing 
ionization equilibrium. The measured log A^(C Ii/C i) is indicated by the red dashed line with the 0.15 dex error in green. Bottom Right: 
The final solutions are shown in black (1 a) and yellow (2 cr). 


Table 4. VPFIT results for the C i data from VLT/UVES 


Comp. 

^abs 

b 

km s“^ 

log N{C i) 
cm“^ 

log N{C I*) 
cm“^ 

log N{C I**) 
cm“^ 

K 

(/1J2) 

1 

2.116069 

7.42±2.35 

12 .88±0.10 

12.66±0.25 

12.54±0.38 

14.6 

(0.29,0.22) 

2 

2.116305 

8.18±0.62 

13.94±0.04 

13.55±0.05 

12.92±0.17 

11.8 

(0.27,0.06) 

3 

2.116425 

4.77±0.99 

13.63±0.07 

13.54±0.05 

12.85±0.16 

18.0 

(0.41,0.09) 

4 

2.116533 

2.62±0.33 

14.55±0.32 

13.90±0.06 

13.61±0.07 

9.1 

(0.17,0.09) 

5 

2.116651 

2.05±0.36 

13.96±0.26 

13.70±0.07 

13.07±0.09 

13.9 

(0.33,0.08) 

6 

2.116808 

7.13±0.92 

13.33±0.06 

13.43±0.05 

13.20±0.07 

27.0 

(0.42,0.25) 

7 

2.117003 

8.05±0.55 

13.43±0.06 

13.80±0.03 

13.51±0.04 

101.4 

(0.52,0.26) 

8 

2.117418 

15.34±2.54 

13.53±0.08 

13.41±0.08 

12.91±0.16 

17.2 

(0.38,0.12) 

9 

2.117625 

5.38±0.47 

13.48±0.07 

13.66±0.04 

13.40±0.04 

35.7 

(0.45,0.26) 

total 



14.84±0.17 

14.56±0.02 

14.18±0.03 

13.5 

(0.30,0.16) 
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factor. We also include the cosmic ray background which 
is known to influence the ionization and chemical state of 
high-density ISM. The cosmic microwave background at the 
absorber’s redshift is also taken into account. The chemi¬ 
cal abundance is initially set to that of the local ISM but 
the gas phase and grain abundances can be increased or de¬ 
creased by a scale factor. We run a grid of simulations for 
different values of hydrogen volume density n(H l) (equiva¬ 
lent to varying the ionization parameter) and subsequently 
vary the above-mentioned parameters to match the observed 
column densities of Zn ll and Fe ll in each corresponding 
velocity component since these two ions have the most reli¬ 
able individual-component measurements from the available 
VLT/UVES data. However, these ions are not sensitive to 
n(H l). We include neutral chlorine and magnesium in the 
modeling of which the column densities are quite sensitive 
to the hydrogen density. The best model and the best-fit 
n(H l) are determined by N{C l) in combination with those 
ions. 

Fig. 1151 shows the Cloudy model for component 3. The 
preferred model is selected not only because the predicted 
N{Zn ll) and A'^(Fe ll) match the measured values but also 
because Mg I, C I, and Cl I agree on n(HI). The models 
that do not satisfy these criteria are therefore rejected. As a 
result, the preferred model gives a metallicity half of that of 
the MW ISM and a radiation field 10 times the intensity of 
the MW ISM radiation. As is often the case, the mo deling 
predicts too much C ll iIFox. Richter fc Fechnerll2014|j . The 
predicted neutral carbon column density increases as n(H l) 
rises. The best-fit n(H l) value can be determined as the in¬ 
tercept of the solid (predicted) and dashed (measured) lines: 
log n(H l) ~ 2. The predicted Cl l does not match the ob¬ 
served column density (Sec. 16.411 until log n(H l) is greater 
than 1.5. Mg l indicates a hydrogen volume density of 
log n(H l) = 2.3. They all agree on a high density (n(H l) ~ 
100 cm“®) neutral medium. The photoionization model with 
this n(H l) value indicates that the electron temperature Te 
~ 70 - 110 K. The average values throughout the slabs are 
consistent with those derived from the C l technique, which 
is expected. 


6.4 Neutral Chlorine and molecular hydrogen 

We detect a strong Cl l A1347 line with seven velocity com¬ 
ponents (Fig. Hill from which we derive a total column den¬ 
sity of log N{Cl l) = 13.95 ± 0.10 cm~^. We disregard the 
region beyond -|- 50 km s“^due to the severe contamination. 
Neutral chlorine and molecular hydrogen are kno wn to have 
a tight relation in the cold phase of t he local ISM (ljur 3 1974 
iMoomev. Federman fc Sheffeij[^012|j . In the presence of H 2 , 
singly ionized chlorine converts rapidly through the exother¬ 
mic reaction Cl^ -I- H 2 —> HCl’^. Neutral chlorine can be 
quickly released from several proces ses especially co l lisions 
between HCl’^ and H 2 or electrons. iBalashev et al.l (l2015l l 
studies this relation at high redshift and finds a ~ 5 (j corre¬ 
lation between A^(C1 l) and A(H 2 ): N{Cl l) « 1.5 x 10“® x 
A^(H 2 ), which confirms that neutral chlorine is an excellent 
tracer of molecule-rich gas. We do not directly detect H 2 
due to a very low S/N but it essentially exists. Using that 
relation, we estimate a molecular hydrogen column density 
of log A^(H 2 ) ~ 19.74 cm“^ and a molecular fraction of 0.1. 



log n(HI/cm 


Figure 15. The Cloudy photoionization model for component 3 
shown in the column densities of atoms/ions vs. n(H l) diagram. 
The measured values are marked with the dashed and dotted lines 
while the predicted values from Cloudy are the solid lines with 
the corresponding colors. The metal abundances are scaled to half 
of the MW ISM. The incident radiation field has an intensity 10 
times that of the MW ISM radiation. 

Higher S/N spectra would be required to show the H 2 ab¬ 
sorption lines. 

6.5 Carbon monoxide 

The VLT/UVES spectrum covers a series of CO A-X bands. 
We have simultaneously fitted 6 CO bands, 0-0, 1-0, 3-0, 4-0, 
5-0, and 6-0, with a two-component (i.e., velocity component 
3 and 4) model. The rotational levels of J = 0, 1, 2 are taken 
into account. As shown in Fig. 1161 individual CO bands are 
very weak but the feature is detected in the stacked spec¬ 
trum. This is the second detection with the simultaneous 
presence of a 2175 A bump beyond the redshift of 2 along 
a quasar line of sight. The total CO column density is log 
N{CO) = 14.1 ± 0.1 cm~^ with an excitation temperature 
of ~ 12 K. fFig. I17II . While the presence of CO is revealed 
thanks to several coincident absorption lines at the expected 
position of the electronic bands, and highlighted by stacking 
the lines, we caution that further observations are required 
to derive the parameters accurately, in particular the exci¬ 
tation temperature. 

The simultaneous presence of CO and the 2175 A 
bump indicates that an efficient shielding of the UV radi¬ 
ation field and high metal abundances may be necessary for 
small carbon-rich grains to exist. These are the ideal condi- 
tions for forming molecules efficiently and in particular CO 
dSloan et al.ll200a ~). 


7 SUMMARY 

We have reported the detection of a strong Milky Way- 
type 2175 A dust absorber at 2 : = 2.1166 towards quasar 
J1211-I-0833 from the BOSS DRIO. Rich metal absorp¬ 
tion lines are detected in the Keck/ESI spectrum and the 
VLT/UVES data reveal a complex velocity structure. It is a 
remarkable system that has simultaneous detections of C l 
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Figure 16. The CO A-X bands with a two-component fit (i.e., 
component 3 and 4). The rotational levels of J = 0,1,2 are in¬ 
cluded in the simultaneous fitting. The last panel shows the 
stacked CO spectrum. 


and CO at the redshift of the 2175 A bump, which are uti¬ 
lized to derive the physical conditions in the absorbing gas. 
Using the Cl fine structure lines and the C ll/C l ratio, we 
provide constraints on the neutral hydrogen volume density, 
temperature, and hence pressure. The results of the analysis 
are: 

1. This 2175 A bump has a bump strength typical 
of a Milky-way-type bump but is wider than the major¬ 
ity of MW bumps. The distribution of MW bump widths 
is non-Gaussian with a mean peak and a tail towards 
larger values, which might s uggest two populations of bumps 
llFitzpatrick fc Massal [200^ 1 . The bump in J1211-1-0833 be¬ 
longs to the population in the tail. If PAHs are the carriers, 
the variati on may be i nterpreted as different mix of PAH 
molecules llDraindl200l . 

2. We measure column densities for a large number of 
metal lines (e.g., Zn ll, Fe ll. Mg ll. Si ll, Ni ll, Mg l, A1 ill) 
covered in Keck/ESI and VLT/UVES. We obtain a hydro¬ 
gen column density of log N{¥L l) = 21.00 ± 0.20 cm“'^ from 
fitting the damped Ly-a profile. With N{}i l), we derive the 



Figure 17. Confidence contours in the vs. log 7V(CO) plane 
assuming 6 = 1 km s“^. The CO column density is log N{CO) 
= 14.1 ± 0.1 cm“'^ with an excitation temperature of ~ 12 K. 


absolute chemical abundances and the dust depletion pat¬ 
tern which suggest that this absorption system has a solar 
Zn abundance and is highly dust depleted with a depletion 
pattern resembling that of the Milky Way disk clouds. 

3. The detailed velocity structure is revealed by the 
high-resolution VLT/UVES data: a minimum of nine ve¬ 
locity components are required to optimally fit the data and 
they span a velocity interval of ~ 150 km s“^. The pro¬ 
files exhibit an edge-leading asymmetry with the strongest 
feature in the blue edge. The asymmetry shown in the pro¬ 
file is consiste nt with the prediction by a simple model of a 
rotating disk (IProchaska fc Wol^ll99^ . 

4. Higher resolution observations with the VLT/UVES 

also reveal the presence of neutral carbon and CO molecules 
in the absorbing gas. The C l fine structure line analysis 
suggests that the physical conditions in the absorbing gas 
are consistent with a canonical (T ~ 100 K) cold neutral 
medium. C l traces dense pockets of cold gas at slightly 
higher pressures. We have further constructed a Cloudy pho¬ 
toionization model which utilizes atoms/ions other than C l 
to constrain the physical conditions in the absorbing cloud. 
Neutral carbon, chlorine, and magnesium mutually agree on 
a hydrogen volume density of n(H l ) ^ 100 cm~^. _ 

5. It has been suggested (e.g.. lElfasdottir et al.l l2009l: 


iNoterdaeme et ^ 12001^ : IProchaska et al. 20091) that a cor¬ 

relation exists between the presence of the 2175 A bump 
feature and prominent C l. The remarkable aspect of this 
quasar absorber is that carbon monoxide is simultaneously 
present along with C l and the 2175 A bump. Since the 
carriers of the 2175 A are generally believed to be carbona¬ 
ceous material that requires neutral carbon and mole cules 
for grains to form and grow (iHenning fc Salamalll99^ . the 
simultaneous presence of C l, CO, and the 2175 A bump 
would not be surprising but rather preferred. Combined with 
the high metallicity, high dust depletion level, and overall 
low ionization state of the gas, all the evidence supports 
the scenario in which the presence of the bump requires an 


















































































J1211 17 


evolved stellar populat ion (i.e., AGB stars) dNoll et al.ll2007l : 
lEli'asdottir et al1[2009l l. The host of the 2175 A bump to¬ 
wards J1211-I-0833 is likely to be an evolved and chemically- 
enriched disk galaxy. 
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